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Abstract 

Time-dependent quantities are calculated in the linear response limit for a correlated one dimen¬ 
sional model atom driven by an external quadrupolar time-dependent field. Besides the analysis of 
the time-evolving energy change in the correlated two-particle system, and orthogonality of initial 
and final states, Mehler’s formula is applied in order to derive a point-wise decomposition of the 
time-dependent one-matrix in terms of time-dependent occupation numbers and time-dependent 
orthonormal, natural orbitals. Based on such exact spectral decomposition on the time domain, 
Renyi’s entropy is also investigated. Considering the structure of the exact time-dependent one- 
matrix, an independent-particle model is defined from it which contains exact information on the 
single-particle probability density and probability current. The resulting noninteracting auxiliary 
state is used to construct an effective potential and discuss its applicability. 

PACS numbers: 31.15.ec, 03.65.-w, 03.67.-a, 71.15.Mb 
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I. INTRODUCTION AND MOTIVATIONS 


There are quantum mechanical problems in which the Hamiltonian depends explicitly on 
time, for example, the interaction of the system with an external time-dependent field where 
H(x,t ) = Hq{x) + H'(x,t). In such cases, the system does not remain in any stationary 
state, and the behavior of it is governed by the time-dependent Schrodinger equation 

d${X,t) 


i h ■ 


dt 


= H(X,t)*(X,t). 


The problem will now be an initial-value (Cauchy) problem of solving this partial differential 
equation for a prescribed initial condition at to- When problems of this short are discussed 
formally, it is common to speak of the perturbation (. H ') as causing transitions between 
eigenstates of Hq. If this statement is interpreted to mean that the state has changed from its 
inital value to a final value, than it is incorrect. The effect of a time-dependent perturbation 
is to produce a nonstationary state, rather than to cause a jump from one stationary state 
to another which were determined by considering boundary conditions solely. 

Experiments open the possibility to investigate dynamical properties of confined systems. 
In that field, of particular interest is the response of the system to time-dependent variations 
of the confining field. Then, even assuming complete isolation of the many-body confined 
system from the environment, there is the question of how the correlation properties change. 
Furthermore, one can consider an atom in its ground state. If at times t > to it becomes 
subject to a time-varying potential energy caused by a charged heavy projectile passing by, 
its electron cloud will be shaken up so that its energy increases. The energy change, related 
to the stopping power, is a measurable quantity. However, when we have shake-up processes, 
the application of one-particle auxiliary pictures is not necessarily useful a priori. 

Motivated by such problems, in this work we consider a simplified interacting model 
system introduced by Heisenberg [l-sj] in the early days of quantum mechanics 


Hq{x!,x 2 ) = ~ 


K 2 ( d 2 


d + TT^oO? + x i) - 7T A ^o0i ~x 2 ) 2 , (1) 


2m \dx 2 dx 2 J ' 2 " uv 1 ' ~ z/ 2 

where A G [0, 0.5] measures the strength of repulsive interparticle interaction. We will use 
ideas and methods 4n6] for solving the quantum motion of a particle in a harmonic oscillator 
with time-dependent frequency, by adding a quadrupolar (Q) time-dependent perturbation, 
of H'(t — y —oo) = H’(t — > oo) = 0 character, to the above ground-state Hamiltonian 


m 


H'(x 1 ,x 2 ,t) = y Q F(t) [x\ + x 2 2 ). 


( 2 ) 
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In the field of heavy-particle interaction with atoms, such term could mimic the shake-up 
process due to close encounters. There, as was stressed [?] by Fermi, quantum mechanics is 
needed since Bohr’s classical treatment is valid [§, 9] only for dipolar perturbation. 

One might think that a two-particle model is a bit trivial to test time-dependent many- 
body methods. But, this is not the case. Indeed, it is the few-body correlated dynamics 
that serve [l(3] as benchmark for methods beyond, for instance, Time Dependent Density 
Functional Theory, TD-DFT. Furthermore, to the theory of breathing modes of many- 
body systems in harmonic confinement, it was utilized 11} to replace a given (say, with 
Coulombic interparticle interaction) Hamiltonian with a quadratic Hamiltonian for which 
an analytic solution exist. A mapping at the Hamiltonian level could be a practical version 


of the isospectral deformation discussed earlier 


12] in the light of the, unfortunately, formal 


character of pure existence theorems and mapping lemmas behind TD-DFT. 

We believe that experience about relevant details in interacting systems is best governed 
by studying time-dependent quantities based on exactly solvable instructive examples. In 
particular, the simultaneous treatment of the notorious kinetic and interparticle energy 
components could be more accurate than in TD-DFT where both are approximated. Besides, 


our result allows a transparent implementation of the mapping-formalism 


13 


14] designed 


to construct independent-particle potentials (and to discuss 15} an alarming paradox arising 
there) by using exact probabilistic quantities of the correlated model as inputs. 

As a last motivation, we note that there have been efforts to clarify the basic relations 
between interparticle correlation and information-theoretic measures for entanglement. For 
instance, the stationary ground-state of the correlated model applied in the present work 
has already been investigated in this respect. A surprising duality between Renyi’s entropies 
characterizing entangled systems, with attractive (A < 0) and repulsive interparticlc interac¬ 


tions was found 


16 


18] and explained 19] recently. The question of how such a remarkable 


duality will change when the correlated system is perturbed by a time-dependent external 
field is an exciting one of broad relevance. One of our goals here is to provide an answer. 

We stress that we restrict ourselves to the linear-response limit, by taking our QF(t ) 
sufficiently small or sudden, since the main goal is on the time-dependent spectral aspects 
of correlation-dynamics and not on comparison with data. However, we should note that 
similarly to the Born series of stationary scattering theory, an order-by-order expansion 
could result in only an asymptotic series. Indeed, higher-order response to external field is 
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an important topic in the energy change during shake-up dynamics of a nucleus 20] or an 


atom 


2l|. The idea that by renormalization of the kinetic energy one could 22] extend the 


validity of linear-response formalism is also a challenging one. Such questions need accurate 
solutions in linear and nonlinear responses. The nonlinear version of the exact determination 
of energy changes in our correlated model system will be published separately. 

The rest of the paper is organized as follows. Section 11 contains our theoretical results. 
Section III is devoted to a short summary and few relevant comments. We will use natural, 
rather than atomic, units in this work, except where the opposite is explicitly stated. 

II. RESULTS AND DISCUSSION 


By introducing the normal coordinates X\ = (op + x 2 )/V2 and X 2 = ( X\ — x 2 )/V% 
can easily rewrite [2] the unperturbed Hamiltonian in the form 


h 2 


H 0 (X 1 ,X 2 ) = - — 


d 2 


+ 


d 2 


Tfl o T 7-0 2 ■\s-2 

+ y -^1 + 7 T w 2 ^-2 5 


2m \dX 2 dX 2 J ' 2 ” 1 “ 1 ' 2 


one 


( 3 ) 


where uj\ = uq and co 2 = (Uo\/l — 2A denote the resulting independent-mode frequencies. It 
is this separated form which shows that a time-dependent perturbation of dipolar character, 
D(t)(x i + x 2 ) = [\/2 D(t)\Xi, will couple only to one normal mode characterized by the 
unperturbed angular frequency cui = u 0 . Therefore, such perturbation jg, 9, 20, 21] produces 
an energy change independent of the correlated aspect of our model. 

In the case investigated here, we have QF{t ) (x 2 + x 2 ) = QF(t)(Xf + X|), i.e., there 
is time-dependent quadrupolar perturbation in both independent normal modes, which 
will evolve in time independently. Thus, we proceed by one oscillator [h(X, f)] in a time- 
dependent harmonic confinement, following the established theoretical path, where one has 
to solve a time-dependent Schrodinger equation, of the form given by Eq. (1), with 

-2 m ?T7 

( 4 ) 




0,1s 


The solution rests on making proper changes of the time and distance scales [4], 5] in order 
to consider frequency variations in fi(f) as it changes from Ho during time-evolution. The 
exact mode, a nonstationary evolving state ^[X, H 0 , i?(h2o), t], contains these scales as 


V>( x,Q 0 ,R,t ) = 


m h2 0 

1/4 

X 2 m Ho 

(, mmX 

hR 2 (t ) 7r 


2 KR 2 {t) 

{ 1 n 0 )_ 




( 5 ) 
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This exact quantum mechanical, time-dependent solution is obtained 4, 5] by considering 
a classical equation of motion 

m + mm = o, (6) 

with complex £(f) = R(t) exp[7 7 (f)] substitution, where the real R(t) gives the length scale 
at the instant t. The nonlinear differential equation, determining this scale becomes 


R + R{t)n 2 (t) = 


(7) 


m)r 

after taking 7 (t) = Q 0 /R 2 (t). The initial conditions are R(t 0 = 0) = 1 and R(t = 0) = 0. 
Since we are dealing with a weak external perturbation in this work on spectral dynamics, 
we solve Eq. (7) via the substitutions R(t) = 1 + r(t) and R{t ) = r[t). 

With r <1 to Eq. (7), we get a forced-oscillator-like differential equation 


f(t ) + (2 = QF{t ), 


( 8 ) 


which is treated by going to the complex notation w{t) = [r(t) + i(2f2 0 ) r (7)] in order to 
derive a first-order differential equation for w(t). Remarkably, Eq. (8) shows transparently 
that the breathing frequency 11] is precisely 2 x Oo, i- e - ; the double of the confinement 
frequency characterizing the ground-state mode of the unperturbed h 0 (X). 

Once an explicit form for F(t) is prescribed to be used in both ( i = 1, 2) differential 
equations for r ? ;(t) [considering the two independent modes, with f2 0 —^ > nq and Oo —^> w 2 ] 
and the solutions for t) are found, the time-dependent wave function becomes 


T(Ai, A^ 2 , t) — cui, R\, t) (U 2 , i? 2 , t), 


(9) 


which is valid for \Q\ < oji, and to which the form for 'ijj(Xi) is given by Eq. (5). By using this 
normal-coordinate-based representation for the exact wave-function we can easily calculate 
with it and the unperturbed Hamiltonian, the expectation values of the kinetic and potential 
energy components, < K(t ) > and < V(t) >, respectively. We get for these quantities 


< K{t) >=Y, 

i=l 


hui 

~4~ 


R?(iUi,t) 


hujj 


+ 


Ri(uJi,t) 


UJ? 


<v(t)>= Y, 


( 10 ) 


(ii) 


2—1 
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In the linear-response limit, the total energy change A E(t) = E(t ) — E(to) becomes 


a m = E =*• E 




i =1 


2—1 


{2ui) 2 +f?(ui,t) 


UJt 


( 12 ) 


To arrive at the consistent r.h.s. above, which is valid for weak external fields, we linearized 
Eqs. (10-11). In this case [with Eq. ( 8 )] one gets Ei(t) = (fauj/4)[2 t) Q F{t)/uf] for 
the rate. Next, we derive for the overlap 0(t) = |'k(t 0 )'I / (f )| 2 the following expression 

1 


0(0 = 


(13) 


\/l + AE 1 (t)/fka 1 y/l + A.E 2 (t) / hu 2 

in terms of AE, (t). This quantity is a measure of correlated (when u\ 7 ^ 002 ) dynamics. Its 
product form reflects the independence of two normal modes during propagation. When 
we are close to the stability limit at A —» 0.5, even a weak external field could result in an 
almost perfect orthogonality at certain, F(f)-dependent, time t > t 0 . 

Now, we turn to quantities which show explicitly the entangled nature of the correlated 
two-particle system in the time domain. By rewriting the wave function in terms of original 
coordinates, we determine the reduced single-particle density matrix [Ti(xi, X 2 , £)] from 


Ti(Xi ,x 2 ,t) = / dx 3 '&*(x 1 ,x 3 ,t)'&(x2,X3,t), 

J —OO 

After a long, but straightforward, calculation we obtain 

T 1 (x 1 ,x 2 ,t) = 0 a (xi ,t) (/>*(x 2 ,t)e~^ A(t)(xi ~ X2)2 , 

where, for further clarihcations below, we introduced the following abbreviations 

mu s (t) 11//4 


(14) 


(15) 


<j>s(x,t) = 


hn 


-S«,(t)* a [l-ia(t)/w,(t)] 


e 2fi 


a(t) 

UJ s (t) 


1 

2 


R 2 (uj 2 ,t)R 2 (uj 2 ,t) 


u 1 


oj 2 


(16) 


(17) 


^ _ 1 [uii(t) — uj 2 (t)] 2 + [Ri(ui, t)/ R\(u\, t) — R 2 (c02 ; t)/ R 2 (u 2 , t)] 2 

4 cui(t) + u 2 (t) 

with c a s {t) = 2(jJi{t)(jj 2 {t)/[ujiit] -E 0^2 (^)], where cUj(t) = tUi/[Ri(uJi, t)] 2 . The one-matrix is 
Hermitian, as it must be, since Ti(xi,X 2 ,t) = ^i(x 2 ,xi,t). Later we will give, as in the 
stationary case [3] earlier, a direct spectral decomposition of this important Hermitian 
matrix in Eq.(15), without considering, as more usual, an eigenvalue problem. 
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The diagonal of the reduced single-particle density matrix gives the basic variable of 

TD-DFT, i.e., the one-particle probability density 

1 


n(x, t ) = 


hn 


e -f 


23 


24] from 


(19) 


( 20 ) 


( 21 ) 


The single-particle probability current, is calculated as usual 

2 hf d 

j(x,t ) — — Re 'S/*(x,x , ,t)——'S/(x,x',t)dx'. 
m J ' i ox 

which, in terms of the functions introduced above, becomes 

... 2 h 

j(x,t) = — xnlXjt) a{t). 
m 

These probability densities satisfy the continuity equation: d t n(x,t ) + d x j(x,t) = 0. 

By taking the A(t) —> 0 substitution in Eq.(15), i.e., neglecting the important role of the 
interparticle coordinate, we can define an auxiliary independent-particle model, where the 
particles move in a certain external field. The wave function of this modeling becomes 

^f s (xi,x 2 ,t) = <j> s (xi,t)(/> s (x 2 ,t). (22) 

This approximate state will still result in the exact probability density and exact probability 
current, by construction. So, one may consider it as an output of TD-DFT. Indeed, using 
Eqs. (6.50-6.52) of [13J] or Eqs. (21-22) of |6], an effective potential V s (x,t) is given by 


23 


/-r d 2 




X 


(23) 


24], In the so-called adiabatic treatment in 


V s [x,t,uj a (t)] = ~uj 2 s (t)x 2 
at least upto a time-dependent constant 

n 

TD-DFT, one uses [13| only the first term as a time-dependent external field. 

But, with such a single-particle potential, which produces from a time-dependent wave 
equation the ^/ s (xi, X 2 , t) state, we also arrive at a dilemma. We know the exact Schrodinger 
Hamiltonian. Especially, we know how the unperturbed Hamiltonian is modified by adding 
a time-dependent quadrupolar perturbation to it. Furthermore, since the perturbation acts 
only over a limited time in our problem, in the Schrodinger picture we can follow the standard 
quantum mechanical prescription to determine the energy change by calculating expectation 
values with the exact wave function T(Xi, X 2 , t) and unperturbed Hamiltonian Hq(Xi, X 2 ). 


In the light of V s \x, t, (u s (f)], the situation in the orbital version of TD-DF 
For instance, it will contain a D s (t)-proportional higher-order derivative 
unfortunately, in TD-DFT, which is based I 3 -I 5 J] on a density-potential mapping lemma, 


1 is not so clear. 
13, 15] . Quite 


we do not 


12] know universal functionals to determine an energy change via them. 
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After the above remarks on open questions in applied DT-DFT, and motivated partly by 


Dreizler’s early suggestion 


12] on using an isospectral deformation (d) and a constrained 


search, we turn to the decomposition of our reduced one-particle density matrix, and, in 
addition, put forward an idea on the possibility of using an approximate nonidempotent one- 
matrix, rf(aq, x 2 , t), designed below. Following our experience Jisj] on the decomposition of 
the stationary Ti^i,^), we will use Mchler’s formula 25, 2b| now to Yi(xi,x 2 ,t)- 

Thus, via the [co s (t) + A(t)] = co(t)[l + z 2 {t)\/[ 1 — z 2 (t )] and A{t) = cu(t)2z(t)/[I — z 2 (t )] 
correspondences, we introduce two new variable uj(t) and z{t). The solutions are 

1 + z(t ) 


U)(t) = 0J 3 (t ) 


1 ~z(ty 


MO + / 

-j 

hf) + 1] 

2 

- 1 

lm \ 

V 

Uw J 



z(t) = 

Now, in term of these we get our point-wise decomposition directly 

OO 

Y 1 (x 1 ,x 2 ,t) = ^2^ ~ z(t)][z(t)] 1 (t)*(x 2 ,t), 

1=0 


(24) 


(25) 


(26) 


where the cpi(x,t) elements of the orthonormal basis set (natural orbitals) are 


M X A) = 


muj(t) 

hro 


1/4 


i 


e 2ft 


V¥T\ 


e 2 R w hA 2 H m {-\/muj{t)/hx). 


(27) 


Of course, extracted quantities, like the exact one-matrix r^aq, x 2 , t) and the associated 
exact probability density n(x,t), contain less and less information then the wave function 
d >(xi,x 2 ,t ) and the pair-function n 2 (xi,x 2 ,t) = ^(xi,x 2 ,t)^f*(xi,x 2 ,t). However, even by 
these extracted quantities one can calculate the exact kinetic, and the external-potential 
energies, respectively, while in TD-DFT only the last energy is, in principle, exact. Thus, to 
formulate an idea, which is based on inversion from known probability density n(x, t ) and 
probability current j(x, t) (i.e., from the crucial fl5 . 27] a(t) phase), we make a rewriting 


n(x,t) = 1 1 ~ z(t)][ z {t)] l \(l>i[xAu(t)\] 2 = ~ z d (t)][z d (t)} l \(pi[x,t,uJd(t)}\ 2 , ( 28 ) 

1=0 1=0 

which is valid if and only if oJdif) = 0J s (t) [1 + Zd(t)\/[ 1 — Zd(t)\. 

So, with one-particle probabilistic functions as input, we can prescribe, from the r.h.s. 

of the above equation, a trace-conserving rf(aq, x 2 , t) which has the form of Ti(xi, x 2 ,t)- 

Such nonidempotent form contains two (interconnected) parameters 0 Jd(t) and Zd(t), instead 






















of the uniquely derived c d(t) and z(t). Whether this deformed one-matrix conld result in, 
because of its tunable flexibility, an essentially better approximation for the kinetic energy 
than the TD-DFT [z(t) = z^{t) = 0] where the one-matrix is idempotent, requires future 
investigations within the framework of a constrained search Jl2, 2s| . 

The information-theoretic aspects of iterparticle interaction, i.e., of a crucial term in the 
Hamiltonian, may be considered via entropic measures which do not use expectation values 
with physical dimensions but use pure probabilities encoded in the occupation numbers 


m = [i -z{t)][z(t)] 1 , 


(29) 


which are, in the present case, time-dependent quantities as well. First, we calculate Renyi’s 


entropy 


29], with q ^ 1 and q > 0, from 


S R (q,A,t) = 


In 


[1 - z(t)\“ 


1-9 1-[«(<)]*' 

The von Neumann entropy is obtained as a limiting case when q —y 1. It is given by 


(30) 


S N (A,t) = — ln[l - z(t)\ - 


z(t ) 


ln[*(*)]- 


(31) 


1 - z(t) 

In the knowledge of z(t), these entropies characterize the deviation from the stationary case 
described by z(t 0 ), i.e., without external perturbation on the correlated two-body system. 
By considering repulsive [A G (0, 0.5)] and attractive [A < 0] versions, we can investigate the 


duality 


18 


19] problem with respect to entropies, now in the time domain. Notice, that the 


knowledge of Renyi entropies for arbitrary values of q provides, mainly from information- 


theoretic point o 


view, more information than the von Neumann entropy alone, since one 


could extract 30] from them the full spectrum of the one-matrix. But, as Eq. (26) shows 


transparently, to physical expectation values one also needs the proper basis set. 

Armed the above, strongly interrelated, theoretical details on the dynamics of interpar¬ 
ticle correlation we come to their implementation. We use a simple model for illustration, 
with H'(t) = Q F{t ) = Q 0(f) exp(— ft1). By tuning the value of the effective rate of change 
(denoted by f3) we can easily investigate the sudden and adiabatic limits. We get 


Ri(A,Ui,t) — 1 + 
-Ri(A, (u^,f) 


Q 


( 3 2 + ( 2^) 2 L 
Q 


e ^ — cos(2o ^t) + sin(2a;j£) 

2a;,- 


P 2 + (2a;;) 2 


[—fie ^ + (2uji) sin(2a;jf) + ft cos(2a;jf)] 


(32) 

(33) 
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AEi(A,uji,t) = 


hoji 


Q 2 


[l + e~ 2pt -2 e~^ cos(2tUjt)] 


(34) 


4 (3 2 + (2uji) 2 ujf 

The sum (i = 1,2) of AEi(A,Ui,t ) measures the total deviation from the unperturbed 
ground-state energy E(t = 0) = (h/2)(uj 1 +CU2). In the kick-limit, where fit — > 0, we have 
AEi(A,Ui,t) (x huji(Q/3t/oJi) 2 /[/3 2 + (2w,) 2 ]. As expected on physical grounds, the change at 
long time, i.e., taking t — > 00 first, becomes a constant which tends to zero in the ft —>• 00 
sudden limit. In the long-time limit the overlap parameter of Eq. (13), which measures the 
time-evolution of orthogonality between exact states, also tends to a constant value. 

Simple independent-particle modelings of the correlated model system could be based 
on effective (e), prefixed external fields {m/2)uj 2 (x\ + x 2 ) before time-dependent perturba¬ 
tion. In such cases one gets A E e (t = oo,tu e ) = 2 x (huj e /4:)(Q/uj e ) 2 /[/3 2 + (2 uj e ) 2 ] f° r the 
total energy change, following the standard quantum mechanical averaging procedure dis¬ 
cussed at Eq. (23). Two options are analyzed here, and a remarkable result is found, as 
Figure 1 signals. Namely, with u> e —* u> s (t = 0), i.e., with the density-optimal modeling, 
the corresponding energy change tends to the exact one at high enough values of since 
2uj~ 1 (t = 0) = {ux} + cu^ 1 ). Somewhat surprisingly, by using the so-called energy-optimal 


Hartree-Fock 


16] modeling, where co e —> a;o\/l — A, one cannot get agreement with the 


exact result at any value of /3. For illustration, a convenient ratio-function is defined as 

AE e (A,p,t) 


K(A,p,t) = 


(35) 


AB(A,/?,«)’ 

and plotted in Figure 1 as a function of the f3 with fixed, density-optimal, c o e — uj s (t — 0). 
We have checked, using both sides of Eq. (12), that with a Q — 0.2 value we are in the 
linear-response limit, where the energy changes are proportional to Q 2 . Notice, that we use 
standard Hartree atomic units, where m — h — 1, to both Figures of this work. 

As we can see from the Figure, far from the sudden limit, i.e., at a small /3 value, the 
density-optimal modeling underestimates the true energy change in the correlated case. We 
speculate that such sensitivity close to the adiabatic limit could be, at least partially, behind 
yet unclarified controversies in the applied TD-DFT 31, 32] to energy transfer by slow ions 
in wide-gap insulators, where A is positive of course. In Section III, we will return to the 
informations given in Figure 1, considering there the energy changes for negative A. 

At the end of this Section, we turn to illustrations of the spectral aspect of dynamic 
interparticle correlation. We apply von Neumann entropy to get a precise insight into the 
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FIG. 1: Ratio of the energy changes lZ(A,t), as a function of t which is measured in atomic units. 
The solid and dashed curves refer to interparticle repulsion (A = 0.4) and attraction (A = —2.0), 
respectively. We used Q = 0.2 and ojq = 1. The top, middle, and bottom panels refer, respectively, 
to the /3 = 0.25, /3 = 0.5, and /3 = 1.0 values. 


time-dependence of an information-theoretic measure. In the light of facts on Hamiltonian- 
based measures, the energy-change and the wave-function overlap at t —> oo, such analysis 
is desirable. Based on it, one can discuss, and maybe understood, challenging interrelations 
between measures which are based on different emphasizes within quantum mechanics. We 
stress here that the exact energy-changes depend, as expected, on the sign of A. 

In Figure 2 we exhibit 5jv(A, t), by taking uj 0 — 1, Q = 0.2, A = 0.40 and A = —2.0. 


These couplings for repulsion and attraction, respective 


at t — 0, in order to shed further light on duality 


16 


y, are fixed to yield equal entropies 


18 


19], now in the time-domain as 


well. Notice, that in the unperturbed case to any allowed [A e (0, 0.5)] repulsive coupling 
there exists a corresponding attractive (A < 0) one for which the calculated entropies are 
equal. For instance, when we are close to the stability limit in the repulsive (r) case, i.e., 
when A r —* 0.5, the corresponding attractive (a) coupling becomes A a —> (l/2)(2A r — 1) _1 . 

The solid and dashed curves refer to repulsion and attraction, as we mentioned. They 
are plotted as a function of the time measured in atomic units. The rate-parameter of the 
time-dependent external perturbation is described by f3 and measured in inverse units of 
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FIG. 2: Information-theoretic, von Neumann, entropies Sjv(A, t), as a function of the time t, which 
is measured in atomic units. The solid and dashed curves refer to a moderate interparticle repulsion 
A = 0.3, and the corresponding (see the text) interparticle attraction A = —0.75, respectively. 
Three values for the rate parameter /3, measured in inverse time-units, are employed to illustration. 
Top panel: (3 = 2. Middle panel: /3 = 5. Bottom panel: (3 = 10. 


time. In order to discuss tendencies, we used three values for this parameter. As expected, 
at very short times the entropies, similarly to the energy changes, grow in time. After about 
1-2 atomic units (i.e., about 25-50 attoseconds) they start to reflect the oscillating charac¬ 
ter builded into the time-dependent occupation numbers Pi (t ), the primary root of which 
is in the time-evolving correlated wave function. The entropies will keep, in absence of an 


environment-made dissipative 


3 


34] coupling, such oscillating behavior without limitation. 


From this point of view, they signal the time-evolvement of the closed interacting system. 
More importantly, one can see from the illustration that the initial dual character of en¬ 
tropies, fixed at t = 0, will disappear due to a new physical variable, i.e., the time. Finally, 
there are reductions in the oscillation-amplitudes by increasing f3. Clearly, at the so-called 
sudden limit (/3 —>■ oo) for a time-dependent perturbation, they will diminish as expected, 
since the system’s reaction-rapidity is limited by its normal mode frequencies. 

In order to justify the above argument more directly, we take now a somewhat less realistic 
modeling by using an abrupt quench of all interactions at to — 0 in the Hamiltonian. In 
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( 36 ) 


that case, instead of Eqs. (33-34), one gets [5j the following exact expressions 

Ri(A,u h t) = \jl + t 2 u 2 

Ri(h,Ui,t) = tU { ■ (37) 

V 1 + t 

There is no change in the system’s kinetic energy since [l/Ri{oJi, t)) 2 + [Ri(u>i, t)/ujj\ 2 = 1. By 
using time-dependent A(t) and a J s (t), which are needed to Eq. (25), we get z(t) = z{t = 0) 
for all t > 0 since [oj s (t)/A(t)\ = 2[2^/ujiuj2 /(cui — 0 J 2)] 2 ■ Thus, there is no entropy- change 
in the total-quench case. This S(t) = constant character could be a useful constraint in 
practical attempts, with a good starting n(x,t), using the interconnected Zd(t) and u>d(t) 
variables of Eq. (28). One can take cOd(t) =>■ u>d/R 2 (uJd,t), where R(uid,t) = \Jl + t 2 u^. 

III. SUMMARY AND COMMENTS 

Based on the time-dependent Schrodinger equation, an exact calculation is performed for 
the time-evolving energy change, and orthogonality of initial and final states, in a correlated 
two-particle model system driven by a weak external field of quadrupolar character. Besides, 
considering the structure of the exact time-dependent one-matrix, an independent-particle 
model is defined from it which contains exact information on the single-particle probability 
density and probability current. The resulting noninteracting auxiliary state is used to 
construct an effective potential and discuss its applicability. We analyzed the energy changes 
in a comparative manner, and pointed out (see, below as well) the limited capability of an 
optimized independent-particle modeling for the underlying time-dependent process. 

A point-wise decomposition of the reduced one-particle density matrix is given in terms 
of time-dependent occupation numbers and time-dependent orthonormal orbitals. Based 
on such exact spectral decomposition on the time domain, an entropic measure of insepa¬ 
rable correlation is also investigated. It is found that the duality behavior characterizing 
the stationary ground-state may disappear depending on the rapidity of the external time- 
dependent perturbation. At realistic rapidity, the undamped time-dependent oscillations, 
found in an entropic measure, reflect the fact that the correlated system evolves in time. 
In this respect, we have a signal on the excited state, similarly as one can see in the time- 
evolving probability density and probability current. 
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Our first comment is based on results obtained for the energy changes, and plotted in 
Figure 1, by considering the sign of A. In the case of an attractive interparticle interaction, 
the case of a nucleus [2], the exact result for AE(A < 0,/3,f) becomes bigger at moderate 
(3 values than the result obtained within the density-optimal framework AE e (A < 0 ,(3,t). 
So, we speculate that a recent result [35] on the excitation of a deformed nucleus within 
TD-DFT could have the same, underestimating character in its quadrupole channel, as the 


atomistic case has. T 


rus, the true enhancement in excitations over the result based on the 


modeling of Teller 36], might also need further theoretical refinement for that channel. 


In the second comment, we focus on the role of the switching rate (3 in the case of energy 
changes. As a preliminary step to future realistic calculations, here we would like to estimate 
a characteristic interaction-time (T oc /3 _1 ) in classical atom-atom collision. To do that, we 
model the three-dimensional screened interaction via a finite-range (r < R ) potential 


V (r) = 




1 - — 
R 


(38) 


Solving the classical problem [37] for the collision time (T) we get 


T(b, v, R) — 


VR 2 - b 2 Z l Z 2 e 2 1 


In 


yfc^R 2 - b 2 + R+ (J Z 1 Z 2 e 2 /2E 1 )' 


(39) 


2A, c 3 / 2 ^/(Z 1 Z 2 e 2 /2E 1 ) 2 + cb 2 

where E\ = M\v 2 /2 is the energy of a heavy charged (Zie) projectile moving with velocity 
v and colliding at impact parameter b with a fixed center of effective charge Z 2 e. The 
dimensionless parameter introduced is c = [1 + (ZiZ 2 e 2 /R)/Ef, i.e., it is related to the 
ratio of potential and kinetic energies. A closer inspection of this compact (and reasonable 
at short-range solid-state conditions) form shows that practically one may use 

R 


T(b, v , R) ~ a(v ) 


E\ 


1 - 


b 2 


(40) 


v E, + Z 1 Z 2 e 2 /R V R 2 

to get a very acceptable estimation as a function of the heavy projectile velocity v, with which 
ct(r>) is decreasing from 4 (at v —> 0) to 2 (at v —> oo). If we take /3(v ) ~ T~ l {y ) and consider 
the energy change determined above, we get a A E(v) \m\~ 2 ~ [a(v)v] 2 character at 
very low velocities. For the opposite, high-velocity, limit A E{y) ~ [f3{v)]~ 2 ~ [cx(u)r^ 2 is 
the scaling, which also looks quite reasonable physically. However, the detailed work on the 
exactly determined non-linear energy changes is left for a dedicated publication. 

There, and also in Bohr’s pioneering modeling [§, 21] with time-dependent dipole fields 


of a passing charge, the time-scale derived above using collision theory, could play a more 
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quantitative role to understand fine details behind experimental predictions 38, |39] on 


energy losses of slow ions in insulators. Furthermore, in the renewed field of time-dependent 
energy losses, a proper time-scale also could help to analyze further a remarkable theoretical 


prediction 


40] on the strong interplay of nuclear and electronic stopping components of the 


observable total energy transfer. Indeed, few-electron shake-up processes can be dominant at 
close encounters in condensed matter. In such dynamical cases, one-electron approximations, 
with double-occupancy for an auxiliary spatial orbital, could result in inaccuracies. 
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